Atom-molecule conversion with particle losses 
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Based on the mean-field approximation and the phase space analysis, we study the dynamics of an 
atom-molecule conversion system subject to particle loss. Starting from the many-body dynamics 
described by a master equation, an effective nonlinear Schrodinger equation is introduced. The 
classical phase space is then specified and classified by fixed points. The boundary, which separate 
different dynamical regimes have been calculated and discussed. The effect of particle loss on the 
conversion efficiency and the self-trapping is explored. 

PACS numbers: 03.65.Bz, 07.60.Ly 



I. INTRODUCTION 

Association of ultracold atoms into molecules is cur- 
rently an active topic in the field of ultracold quantum 
physics, it attracts much attention due to its important 
applications ranging from the production of molecular 
Bose-Einstein condensates to the search for the perma- 
nent electric dipole moment, see for example jll-ll2l|. By 
applying a time varying magnetic field in the vicinity of 
Feshbach resonance, a pair of atoms can bound into a di- 
atomic molecule [HI, EH, this conversion can be described 
by the Gross-Pitaevski (GP) equations within the mean- 
field theory (MFT) |15l420| . Such an treatment reduces 
the full many-body problem into a set of coupled non- 
linear Schrodinger equations and the complicated many- 
body dynamics is then turned into a two-mode dynam- 
ics. Earlier study shows that the nonlinearity, which 
arises from both the atom-atom and molecule-molecule 
couplings, plays an important role in the system. Four 
distinct regimes, each has different feature in dynamics 
can be classified, accordingly the bifurcation of the fixed 
points in the classical phase space [H, Hoj] is identified. 

Decoherence that arises from the unavailable coupling 
between the environment (thermal atoms or molecules) 
and the condensed system plays an important role in 
atomic or molecular Bose-Einstein condensate [2l| - |26| . 
Description of decoherence by fully including the quan- 
tum effects requires sophisticated theoretical studies, 
however the standard approach in quantum optics can 
reduce the complexity and in fact it has been widely 
used in Bose-Einstein condensates [27- 32] . For an atom- 
molecule conversion system, we then ask: will the deco- 
herence effect be different from that in the pure atomic 
or molecular Bose-Einstein condensates? What are the 
fixed points in this situation? How do these fixed points 
behave? We will answer these questions in this paper. 

In this paper, we focus on the amplitude decoherence 
(particle loss), which may arise from inelastic collision 
between condensate and noncondensate particles in the 
system. The standard approach in quantum optics for 
open systems are used and the master equation is derived 
by treating the noncondensate atoms and molecules as a 
Markovian reservoir. Under the mean-field approxima- 



tion, an effective non-Hcrmitian Gross-Pitaevskii equa- 
tion is derived. Bifurcation of the fixed points divides 
the parameter space into different dynamical regimes, 
boundaries that separate these regimes are changed by 
the decoherence. By calculating the corresponding Ja- 
cobian matrix, we find that a sudden transition in the 
fixed point from elliptic point to attractor or repeller 
happens with non-zero decoherence rate, which reflects 
the meta-stable features of the system under the deco- 
herence. The atom-molecule conversion efficiency for the 
molecular condensate as well as the self-trapping are also 
studied. 

The paper is organized as follows. In Sec. UH we in- 
troduce the model and transform the master equation to 
a non-Hermitian nonlinear Schrodinger equation. Condi- 
tions that determine the fixed points are derived. In Sec. 
IIII1 we define different regimes by the bifurcation of the 
fixed points and study the dynamics of system in these 
regimes. In Sec. IIV1 we investigate the effect of particle 
loss on the conversion efficiency. In Sec. [V] we shed light 
on the self-trapping taking the decoherence into account, 
an explanation for the observed features is given in the 
framework of mean-field theory. Finally, we conclude our 
results in Sec. [VI] 



II. MODEL 

Based on the two-mode approximation, the Hamilto- 
nian that includes the atom- atom collision U aai atom- 
molecule conversion with rate V, and molecule-molecule 
couplings Ubb takes the following form [l6l []~7jj 

H = fjtaO/ a + Hbb'b + U aa a'aJaa + Ubbb^b^bb (1) 
+U ab a*aPb + V^a^b + tfaa). 

The master equation [33j that takes the particle loss into 
account can be derived as in the textbook [2lj . 
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p = — i[H,p] ■^-{a)ap + pa) a — 2apa)) (2) 
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where T a and represent decoherence rates for atomic 
and molecular modes, respectively. In the mean-field ap- 
proximation, the quantum fluctuation is negligible. It is 
appropriate to replace a and b with c numbers a — \a\e a 
and b = \b\e l9b . With these considerations, the master 
equation (|2|) can be casted into the following nonlinear 
Schrodinger equation, 



H = 



R-Uz- 
Va 



;T a 2Va* 
-2R + 2Uz 



(3) 



(4) 



with z = | a. | 2 — 2|6| 2 describing the number difference for 
atoms in the two modes. U — jU ab — \Uaa — \U bb rep- 
resents a coupling strength and V is the conversion rate. 
R = \{2n a - fi b + 2U aa - \U bb ) denotes the energy dif- 
ference between the two modes, which can be effectively 
tuned by a time- varying external field [lj], El] ■ Here and 
hereafter, we rescale R, U, Y a , Y b in units of V, and t in 
units of 1/V, h = 1 has been set, hence all parameters in 
this paper are of dimensionless. 

Similar to the two-mode Bose-Hubbard model, the pro- 
jective Hilbert space for such an atom-diatomic molecule 
conversion system can be spanned by a set of Bloch vec- 
tors. Under the mean-field approximation, the Bloch vec- 
tors can be written as [19( 



h = (2v / 2i?e[(a*) 2 6],2v / 2/m[(( 



*\2 



61, lal- 



2|6| 2 ). (5) 



With the normalization condition \a\ 2 + 2\b\ 2 = 1, the 
projective space is a tear-drop shaped surface as shown 
in Fig. ® 

To analyze the dissipative dynamics of the system in its 
classical phase space, we define relative phase 9, particle 
number n and normalized population difference S as 



S 



= 26 a - , 

= 2|6| 2 4 

z 



(6) 
(7) 

(8) 



It is found that the bifurcation of the fixed points falls 
into four regimes in the parameter space (l6l [20| as shown 
in Fig. [TJ Due to inelastic collisions between particles in 
the condensates and that in the thermal cloud, the loss of 
atoms and molecules as one source of decoherence is un- 
avoidable in practices. Thus it is interesting and desired 
to study the dynamics of this system with decoherence. 
By the standard procedure used in non-Hermitian sys- 
tem [13, HU , we study the dynamics of the system with 
notations C — Un, = Vy/n. With these knowledge, 
the effects of the particle loss on the dynamics as well 
as the features of fixed point can be clearly revealed in 
classical phase space. For simplicity, in the following dis- 
cussion, we set r = r_ and rewrite the equations for the 
population difference and relative phase as 



S = -2Q(1 + S)VT~S sind -T{l-S 2 ), (13) 

(14) 



1 _ Q C 

4CS -AR- n , cos 9. 



VT^S 

The fixed points of the system are determined by 



S = = O. 



(15) 



To be specific, the fixed points on the boundary are S 

— 1, 9 = arccos(— - 
are determined by 



1, 9 = arccos(— ^^jp-^ ), while the other fixed points 



(9r 2 + 64C 2 )5 3 - (r 2 - An 2 + 64i? 2 ) 

-(15r 2 - 36ft 2 + 64C 2 + 128CR)S 2 

-(24ft 2 -7r 2 -64i? 2 - 128CR)S = 0, (16) 
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Through a Jacobian matrix defined by 
( dS/dS dS/89 \ 

V d <i/ ds dd l dd ) ' 



(17) 
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Inserting these definitions into Eq. ([3]) , a set of evolution 
equations is obtained 



S = -2VV^(1 + S)VT^S sin8 ~T-(l- S 2 ), (9) 



= AUnS - AR - V 
= -(T++T-S)n, 
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cos 9, 



(10) 
(11) 



with r + = i(r a 



T b ) and r_ = ^(r a — r^) represent- 
ing the total and relative decoherence rates for the two 
modes, respectively. Ignoring the coupling between the 
system and its environment (r a = T b = 0), the dynamics 
of the system can be described by an effective Hamilto- 



man as 



H = 2V{1 + S)VT^S cos 9- 2US 2 + ARS. (12) 



we can study the stability of the fixed points as did in 
the literature (Hill Hi]. 

From Eqs. (fT3"|) and (fl"4"l). we find that the dynamics 
of the system depends only on the relative decoherence 
rate T, but not on the total decoherence rate. The total 
decoherence rate affect the system through n = noe~ r+t . 
With zero relative decoherence rate (r = 0), the effective 
Hamiltonian can be written as 



H = 2Q(1 + S)VT^S cos 9 -2CS 2 +ARS 7 (19) 

which differs from the case without decoherence only at 
the conversion amplitude ft and coupling strength C. In 
the next section, we shall discuss the effect of relative 
decoherence rate on the dynamics of the system. 
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FIG. 1: (color online) Parameter space spanned by nonlinear- 
ity C and energy difference R. Different regimes are distin- 
guished by boundary lines, where black solid lines represent 
the case for F = 0, while blue dashed lines denotes the case for 
r = 1.2 in (a). Green solid lines in (b) denotes the boundary 
lines for Y = 2.4. (c), (d), (e), and (f) describe the classical 
phase space for region /, //, III, and IV, respectively. 




FIG. 2: (color online) Mean-field dynamics on the Bloch 
sphere for Hermitian (left) and non-Hermitian (right) cases. 
The north pole and south pole of the sphere corresponds to 
the pure atomic condensate and the pure molecular conden- 
sate, respectively. Red spots and center of the vortex denote 
the location of the fixed points. Blue solid lines repersent the 
trajectories for the time evolution of the system. Parameters 
chosen are R = 1, U = for (a) and (b), R — 0, U = 2 for (c) 
and (d), and R = 0, U = for (e) and (f). The spheres on 
left side ((a),(c),(e)) describe decoherence free case (F = 0), 
while the spheres on right side depict the decoherence case 
(T=l). 




FIG. 3: (color online) (a) Locations of fixed points for differ- 
ent decoherence rate F and energy difference R. Parameters 
chosen are U = 0, V = 1, and F = 0,0.9,1.6 for red solid 
line, green dashed line and blue dotted line, respectively, (b) 
Locations of fixed points for different decoherence rate V and 
interaction strength C. Parameters chosen are R = 0, V = 1, 
and r = 0, 0.5, 1.5 for red solid line, green dashed line and 
blue dotted line, respectively. 



III. FOUR DYNAMICAL REGIMES WITH 
DECOHERENCE 

In Ref. [16|, without considering the decoherence ef- 
fect, by the feature of fixed points, the parameter space 
was divided into four regions. Here we re-divide the re- 
gion by taking the decoherence into account (see Fig. Q] 
and Fig. [5]). Boundaries that separate different regions, 
are determined by numerically solving Eqs. (|16I17[) . Note 
that the fixed points on the boundary behave like the 
fixed points in the region labeled by a smaller number 
(e.g. boundary that separate regions I and II belongs to 
the region 7). 

Figure QJc) shows Poincare section of the classical 
Hamiltonian for the region I. The only fixed point is lo- 
cated near the boundary of the phase space (5 = 1) and 
the dynamics of the system is localized. When taking 
the decoherence into consideration, the fixed point near 
5 = 1 turns into an attractor, where Figs. [SJa) and^b) 
show trajectories on the tear-drop shaped Bloch sphere. 
The dynamics of the system becomes delocalized due to 
the appearance of such an attractor. 

By changing the energy difference R and the nonlin- 
earity C (see Fig. Ufa)), the system can go across the 
boundary into region II, the fixed point in region / bifur- 
cates into two elliptic points and a hyperbolic one as Fig. 
[TJd) shows. The region II shares similar features with 
the self-trapping in the two-mode Bose-Hubburd model 
[3IH [37( 1 . With a negative decoherence rate, both of the 
two elliptic fixed points transit to attractors in this re- 
gion (see Figs. [2jc) and[2jd)). While the locations of 
the stable attractors are just slightly changed due to the 
decoherence (see Fig. |3Kb)). 

Figure [He) illustrates the Poincare section of the clas- 
sical Hamiltonian for region III without decoherence. In 
this region, large amplitude oscillations around the ellip- 
tic fixed can be observed, see Fig. [He). With U = and 
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R = 0, the location of the fixed points in this region can 
be derived analytically 




arcsin(-^- 



) 



2vr-arcsin(-^-) 



(20) 



where we assume the decoherence rate positive, and the 
relative phase was restricted in 9 € [0, 2w]. From Eq. 
(j20|) . we find that in addition to the feature changes of 
the fixed points, the relative phase between the two fixed 
points decreases and the fixed points becomes asymmet- 
ric due to the decoherence as shown in Fig. [2Jf ) . As 
the decoherence rate increases, the area of regime III 
is compressed (see blue dashed line in Fig. H{a)). The 
two boundaries coincides and region III vanishes (see 
dash dotted line in Fig. G^a)), when decoherence rate 
is larger than a threshold (r > v2fi)j a hyperbolic fixed 
point arises from the bottom of the phase space (see dash- 
dotted line in Fig. [3](a) ) . The boundary that separates 
regions III and IV is shifted due to decoherence. This 
boundary shift can be explained as a threshold decrease 
in the energy difference R (denoted by Rq and R\ in Fig. 
[3Ja)), which is an witness for the bifurcation of fixed 
points in classical phase space. 

The dynamics in region IV behaves similarly as that 
in region /. The elliptic fixed point turns into an attrac- 
tor due to negative decoherence rate, the dynamics in 
this region then becomes delocalized (see Figs. Ufa) and 
lib))- 

Next, we focus on the changes of the fixed points. Such 
a change in classical phase space is fundamental for non- 
hermitian Bose-Hubbard system [2?], Hi], [HJ . However, 
we find that, in the atom-molecule conversion system, 
the change differs from Bose-Hubbard model in two re- 
spects. Firstly, the type of the fixed point (a repeller or 
an attractor) is determined by the the sign of decoher- 
ence rate F and the location of the fixed point S. If T 
and S are different in sign, i.e., one of them is positive 
while another is negative, the original elliptic fixed point 
transits into a stable attractor. Otherwise, the original 
fixed point turns into an unstable repeller. Secondly, the 
transition is sudden, in other words, the transition hap- 
pens provided the decoherence rate is not zero. This 
is different from the decoherence effect on Bose-Einstein 
condensates in a double-well potential, namely there ex- 
ists a critical value for the decoherence rate [27|. In the 
atom-molecule conversion system, the transition happens 
once the decoherence exists, regardless of how small the 
decoherence being. The phenomenon reflects not only the 
meta-stable behavior of the open many-particle system, 
but also the sensitivity of the atom-molecule conversion 
system to particle loss. 



IV. CONVERSION EFFICIENCY FOR 
MOLECULAR CONDENSATE 

In experiments, the association of ultracold atoms into 
diatomic molecules can be achieved by applying a time- 
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FIG. 4: (color online) Conversion efficiency W in (a) and 
relative efficiency M in (b) as a function of the sweeping rate 
ft under different decoherence rate Y. 



dependent magnetic field in the vicinity of a Feshbach 
resonance, which corresponds to the change between dif- 
ferent regimes (I — > III — > IV) in the parameter space 
(see Fig. [T|). To examine the effect of decoherence on 
the conversion process, we define conversion efficiency, 
relative efficiency and sweeping rate of the external field 
as 



W 
M 



\KT)\ 2 
n{T) ' 

W(T) - W(0) 



W(0) 



= R, 



(21) 

(22) 
(23) 



where T denotes the terminal time for the conversion, 
VF(r) and W(0) denote the conversion efficiency with 
and without decoherence, respectively. M describes the 
relative increases or decreases of the efficiency between 
decoherence and decoherence-free case. By adjusting 
the external magnetic field [13], R can be linearly ma- 
nipulated to across the Feshbach resonance point (R = 
Ro+l3t 7 R — (3T, t £ [0, T]), until the system relaxes into 
a steady state. The conversion efficiency with decoher- 
ence can be calculated with the same parameters except 
the decoherence rate. Here we choose the initial state of 
the system to be in pure atomic mode (|a(0)| 2 = 1). 

The results of W show that conversion efficiency in- 
creases with positive decoherence rate. While a negative 
decoherence rate decreases the conversion efficiency (see 
Fig. 13]). This can be interpreted by the appearance of 
attractor or repeller in the phase space. With a negative 
decoherence rate, the elliptic fixed points near the atomic 
mode turns into an attractor and the atoms are attracted 
to stay away from molecular mode, (see Figs, [^b) and 
HJf)). The conversion process is depressed by such an at- 
tractor and the conversion efficiency decreases. Similarly, 
a positive decoherence rate will increase the conversion 
efficiency. 
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FIG. 5: (color online) Time evolution for the population of 
atomic mode P(a) = |a(i)| 2 under different decoherence rates 
as T = 0, —0.5, 0.5 denoted by red solid line, black dashed line 
and blue dash-dotted line both in (a) and (b). Parameters 
chosen are V = 1, R = for both (a) and (b), (7 = for 
(a) and U = 1.5 for (b). The initial population for atoms are 
|o(0)| 2 = 0.9. 

V. TUNNELING AND SELF-TRAPPING 

In this section, we investigate the effect of particle loss 
on the dynamics of the system, the atoms may oscil- 
lates between atomic and molecular modes (correspond- 
ing to regime III) , and they can also be trapped in one of 
the modes (corresponding to the regime II in parameter 
space). 

In regime III, the atoms oscillate between atomic mode 
and molecular mode (see Fig. [He)). When the relative 
decoherence rate is positive, the fixed point transits from 
elliptic to a repeller, the amplitude of the oscillation is 
then increased (see dash dotted line in Fig. EJa)). While 
for negative relative decoherence rate, the oscillation is 
compressed, since the elliptic fixed point suddenly tran- 
sits to an attractor (see dashed line in Fig. [S^a)). 

With C increases, the dynamics of the system turns 
into the self-trapping regime, which belongs to the regime 



II in Fig. [IJa). We find that the threshold of the coupling 
constant is decreased by the decoherence, i.e., the deco- 
herence supports the self-trapping (denoted by Co and 
C\ in Fig. GUb)). With negative relative decoherence 
rate, the fixed point near the atomic mode transits into 
an attractor. The self-trapping in atomic mode keeps(see 
black dashed line in Fig. [Sfb)). When the relative de- 
coherence rate is positive, which indicates a repeller in 
the phase space, the self-trapping in atomic mode is ru- 
ined, because the atoms are repelled and converted into 
molecules, as dash dotted line shows in Fig. [5jb). 
A physics understanding is coming 

VI. CONCLUSION 

In summary, we have investigated the effect of parti- 
cle loss on the dynamics of the atom-molecule conversion 
system. Within the mean-field approximation, the clas- 
sical phase space is specified and the fixed points are 
calculated. Due to the bifurcation of the fixe points in 
the phase space, the parameter space can be divided into 
different regimes. We find that the boundary, which sep- 
arates different regimes are changed by the decoherence. 
A sudden transition for the fixed points from elliptic 
to attractor or repeller happens. Such a transition not 
only reflects the meta-stable behavior of the system, but 
also characterizes the phase-space structure of the atom- 
molecule conversion system. The effect of decoherence on 
the conversion efficiency and the self-trapping is explored 
with the mean-field approximation. 
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China Normal University, as well as the National Re- 
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